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INTRODUCTION 


The entry trajectory of the shuttle orbiter is designed for rapid deceleration of 
the vehicle at high altitude. The resulting combination of low velocity at high altitude leads 
to significant departures from chemical equilibrium in the flow field for that portion of the 
flight where real gas effects are important. In this paper, the extent, and some of the 
consequences, of nonequilibrium effects in the inviscid orbiter flow field are examined. 

To this end, calculations have been made for the reacting flow along selected streamlines 
in the windward plane of symmetry of a blunted delta shape representative of a shuttle 
orbiter. 

An approximate model of the pressure field in the windward plane of symmetry of a 

blunted delta wing was constructed from wind tunnel measurements of the body surface 

pressure distribution and approximate calculations of the shock shape and streamline 

locations. Then numerical solutions for the flow along streamlines were obtained for 

this pressure distribution. Solutions were obtained for both equilibrium and finite- reaction- 

rate, nonequilibrium chemical reactions. Since the streamtube solutions are much easier 

to compute than fully three-dimensional, reacting flows, results were obtained for a wide 

range of conditions. The validity of this approach is based on the assumptions that the 

pressure field is insensitive to nonequilibrium effects and that the flow along streamlines 

is quasi-one dimensional. Numerical solutions have been used to verify this approach for 
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other blunted slender shapes such as spherically blunted cones. 



The experimental surface pressure distributions were obtained in the NASA Langley 
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Research Center Mach 8 Variable-Density Tunnel. Pressure data were obtained for two 
shuttle test shapes, a blunted delta and a straight body, at angles of attack of 0, 20°, 40°, 
and 60°. All of these experimental results are presented even though the reacting flow 
calculations were confined to angles of attack of 20° and 40°. The streamline flow calcula- 
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tions were obtained with a computer program previously developed at the Cornell Aero- 
nautical Laboratory for the quasi-one-dimensional, reacting flow of a general gas mixture. 

* 

In the present work, air was assumed to be composed of a neutral species N 2 » O^, NO, N, 

O and Ar. Ionization of NO was also included. The thermodynamic and chemical kinetic 
data employed in the calculations are listed in Ref. 3. 

The present paper first presents the data for the body surface pressure distributions, 
the construction of the approximate model of the pressure field in the windward plane of 
symmetry is described. Then results of both equilibrium and nonequilibrium flow calcula- 
tions are discussed for three typical shuttle trajectory points in the altitude range of 60.96 to 
76.20 km (200,000 to 250,000 ft). Calculations are compared for two angles of attack, 20° 
and 40°, and two nose radii, 0.305 m (1 ft) and 1.22 m (4 ft). Finally, the implications of 
the departures from chemical equilibrium to shuttle technology are discussed. The questions 
of surface heating rates, boundary layer transition, and entropy layer swallowing are treated. 

' The nonequilibrium flow studies were supported by the NASA Langley Research Center 
under Contract No. NAS 1-10710. 
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PRESSURE TEST MODELS 
(Fig. 1) 


A series of experiments to measure the windward surface pressure distributions on two shuttle 
orbiter test shapes were conducted in the NASA Langley Mach 8 variable-density hypersonic tunnel, 
face pressures were obtained for the blunted delta and straight bodies shown below. 
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DELTA-WING AND STRAIGHT-BODY MODELS FOR 
WIND TUNNEL PRESSURE MEASUREMENTS 
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TEST CONDITIONS FOR PRESSURE MEASUREMENTS 
(Fig. 2) 


The tests were conducted at a nominal test-section Mach number of 8.0, unit Reynolds numbers in the 
range 8.8 x 10^ m ^ to 41.5 x 10^ m \ and angles of attack of 0°, 20°, 40°, and 60°. The Reynolds num- 
bers given in the figure are based on model length. 
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EXPERIMENTAL TEST CONDITIONS 
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Figure 2 













PRESSURE DATA FOR STRAIGHT BODY 
(Fig. 3) 


In this figure, the data from the experiments on the straight-wing body are 
tabulated. The data are given in the form where .fi is the model 

stagnation pressure. The streamwise distance, S '/ R N , is measured from the model 

centerline, b is the local span. 


CENTERLINE PRESSURE DATA FOR THE STRAIGHT BODY 


RUN NO. 

28 

29 

30 

31 
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35 

36 

s/r n 

P/P T 

'2 

P/P T 

\2 

P/P T 

'2 
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.1331 
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.4983 
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Figure 3(a) 




SPANWISE PRESSURE DATA FOR THE STRAIGHT BODY 


RUN NO . 

28 
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.1594 
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.8518 
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PRESSURE DATA FOR THE DELTA BODY 
(Fig. 4) 

In this figure, the data are tabulated for the experiments on the blunted delta wing. Again the 
data are normalized by the model stagnation pressure and the streamwise distance is measured from the 
model centerline. 


CENTERLINE PRESSURE DATA FOR THE DELTA BODY 
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Z071 
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.1459 
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.4691 
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.8466 

.8314 
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.1184 

.1347 

.4963 
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.4849 
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.1684 
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Figure 4(a) 



SPANWISE PRESSURE DATA FOR THE DELTA BODY 


RUN NO . 

3 

4 

5 

6 

10 

11 
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S / R n 
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p / p t 2 

P / P T 

*2 

p / p t 2 
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2 
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.1618 

.4849 

.4832 

.8339 

.8298 

— 
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-.193 

o 

.1360 

.1516 

.4868 

.4748 

.8514 

.8466 

.8212 


.193 

.1360 

.1488 

.4887 

.4748 

.8514 

.8410 

.8212 


.484 

.1418 

.1572 

.4792 

.4748 

.8299 

.8241 

.8008 

13.61 

-.664 

— 

— 

_ 

— 

— 


— 


-.483 

.1554 

.1684 

.4868 

.4944 

.8397 

.8438 

.8161 


-.302 

.1477 

.1628 

.4887 

.4944 

.8494 

.8494 

.8212 


-.121 

.1516 

.1657 

.4887 

.4973 

.8534 

.8607 

.8264 


0 

.1516 

.1657 

.4887 

.4973 

.8553 

.8579 

.8264 


.121 

.1558 

.1657 

.4906 

.4973 

.8553 

.8607 

.8264 


.302 

.1516 

.1657 

.4906 

.4973 

.8494 

.8523 

.8264 


.483 

.1554 

.1657 

.4906 

.4917 

.8378 

.8410 

.8110 


.664 

.1593 

.1684 

.4811 

.4860 

.8124 

.8184 

.7854 

20.61 

-.798 

.1574 

.1741 

.4697 

.4804 

.7928 

.7847 

.7651 


-.588 

.1516 

.1684 

.4868 

.4860 

.8553 

.8579 

.8212 


-.378 

.1477 

.1628 

.4982 

.4944 

.8612 

.8720 

.8366 


-.210 

.1477 

.1600 

.5001 

.4973 

.8593 

.8720 

.8366 


-.084 

.1457 

.1572 

.5020 

.5001 

.8593 

.8720 

.8366 


0 

.1457 

.1572 

.5020 

.5001 

.8593 

.8720 

.8417 


.084 

.1457 

.1572 

.5020 

.5001 

.8593 

.8748 

.8417 


.210 

.1477 

.1600 

.4982 

.5001 

.8612 

.8748 

.8468 


.378 

.1477 

.1684 

.4924 

.4973 

.8632 

.8748 

.8468 


.588 

.1535 

.1741 

.4849 

.4973 

.8534 

.8466 

.8314 


.798 

.1574 

.1769 

.4678 

.4832 

.7811 

.7819 

.7549 

3 Z 61 

-.856 

.1328 

.1741 

.4639 

.4776 

.7440 

.7481 

.7344 


-.635 

.1516 

.1657 

.4792 

.4944 

.8045 

.8101 

.7906 


-.414 

.1457 

.1628 

.4899 

.5086 

.8202 

.8241 

.8110 


-.249 

.1457 

.1628 

.4811 

.5029 

.8221 

.8269 

.8161 


-.138 

.1457 

.1628 

.4811 

.4944 

.8241 

.8298 

.8161 


-.055 

.1418 

.1572 

.4754 

.4860 

.8241 

.8258 

.8161 


0 

.1418 

.1572 

.4754 

.4860 

.8241 

.8298 

.8161 


.055 

.1438 

.1572 

.4734 

.4889 

.8241 

.8298 

.8161 


.138 

.1399 

.1572 

.4716 

.4917 

.8221 

.8269 

.8110 


.249 

.1438 

.1600 

.4754 

.5001 

.8221 

.8269 

.8110 


.414 

.1438 

.1600 

.4754 

.5001 

.8182 

.8241 

.8059 


.635 

.1477 

.1600 

.4716 

.4832 

.8026 

.8101 

.7957 


.856 

.1554 

.1713 

.4527 

.4720 

.7381 

.7424 

.7344 


Figure 4(b) 


















SCHLIEREN PHOTOGRAPHS OF MODEL FLOW FIELDS 

(Fig. 5) 

In the series of experiments to measure surface pressure on the shuttle test shapes, 
Schlieren photographs were also taken. The photographs from two of the test conditions are 
shown in this figure. The cases illustrated are the blunted-delta and straight body at an 
angle of attack of 40°. 

For the blunted delta wing the shock shapes over the windward surface are quite 
straight for oC = 20° and oC = 40°. In the nonequilibrium flow calculations described later 
in the paper the flow downstream of the nose is assumed to be a uniform flow region and this 
assumption appears to be accurate. 

The comparison of the two Schlieren photographs shown here indicates that the 
difference in the body shapes is reflected in the shock shape for this angle of attack. Both 
the Schlieren photographs and the pressure data indicate that the cross flow begins to affect 
the inviscid flow in the windward plane of symmetry at oc = 40°. At cx = 60°, the 
difference in the shock shapes is more pronounced. Also, the windward centerline pressure 
decays toward the base of the body; the pressure decreases more rapidly for the straight-body. 
At an angle of attack of 20° the shock shapes and windward centerline pressure distributions 
are nearly the same for these two body shapes. 






WINDWARD SIDE CENTERLINE PRESSURE DISTRIBUTIONS 
FOR BLUNTED -DELTA SHAPE 
(Figs. 6 and 7) 

In a previous study'*' of the hypersonic flow over blunted slender cones the surface 
pressure distribution was found to be insensitive to nonequilibrium effects. Since the flow 
along the windward plane of symmetry of the blunted-delta wing should be much like a cone 
flow this result has been used in the present study to scale the measured pressure distri- 
butions to flight conditions. The data were normalized by the tunnel-model stagnation pressure 
and the pressures on the flight vehicle were obtained by using the appropriate value of the 
stagnation pressure. The initial portion of the pressure distribution, i. e. between the 
stagnation point and the first data point, was taken from a composite correlation'* of ideal 

£ 

od gas, equilibrium and nonequilibrium flow calculations, and experimental data for the flow 

over spherically blunted bodies. The stagnation point was assumed to be the point where 
the body surface is normal to the free stream direction. A full scale body length of 
45. 6 m (150 ft) was used. 

The body surface pressure distributions which were used in the streamtube 
calculations for oc = 20° and oc = 40° are shown in the next two figures, together with 
the experimental data on which they are based. Here and throughout the remainder of the 
paper denotes distance along the surface streamline, measured from the stagnation 


point. 


PRESSURE DISTRIBUTION ALONG WINDWARD CENTERLINE FOR 
BLUNTED-DELTA BODY, oC= 20° 



Figure 6 
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PRESSURE DISTRIBUTION ALONG WINDWARD CENTERLINE FOR 
BLUNTED-DELTA BODY, oc= 40° 


RUN 5 
RUN 6 


ASSUMED PRESSURE DISTRIBUTION 
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APPROXIMATE CALCULATION FOR SHOCK SHAPES AND STREAMLINE LOCATIONS 

(Figs. 8 and 9) 

In order to obtain the pressure histories along streamlines off the body surface, an 
approximate model of the pressure field was devised. The shock shapes and streamline 
locations were obtained, for equilibrium flow and again an appeal was made to the results for 
the flow about blunted cones. These results indicate that the shock shape as well as the 
surface pressure is insensitive to nonequilibrium effects. 

The bow shock wave was assumed to be circular in the nose region and straight on 
the afterbody. The final shock angle was determined from the asymptotic, experimental 
surface pressure and oblique shock solutions for equilibrium flow. The circular portion of 
the shock was not taken to be concentric with the body. Rather, the center and radius of 
curvature were found by using a correlation for shock standoff distance at the stagnation 
point^ and a mass balance to get the shock standoff distance at the body shoulder. The mass 
balance was done by assuming the flow in the stagnation region of the windward plane of 
symmetry to be the same as that in the axisymmetric flow over a blunted cone. The circular 
portion of the shock was joined to the straight portion at the point of tangency. In the nose 
region, the streamline locations also were found using a mass balance technique. The mass 
flux was assumed to vary linearly from the shock to the body along normals to the body surface. 

The solution for the equilibrium flow along the body streamline was computed using the known 
pressure distribution and an isentropic expansion from the known stagnation conditions. The 
values at the shock were obtained from the shock angle and shock tables for equilibrium flow. 5 
In the afterbody region the streamlines were assumed to be parallel to the body. This assumption 
of a uniform pressure field on the afterbody appears justified both by Schlieren photographs 
obtained in this study (see Fig. 5) and three-dimensional flow calculations for equilibrium flow 
over a blunted delta wing. Having established the streamline locations, the pressure distributions 
along them were computed assuming a linear variation between the pressure at the body and the 
shock, along normals to the body. 

The shock shapes and streamline locations obtained for c>d = 20° and ot = 40° are 
illustrated in the accompanying figures. The flow field consists basically of two regions: 
those streamlines which are processed by the strong curved shock and those which pass through 
the limiting oblique shock. The streamline which passes through the point where the shock 
becomes straight is referred to as the outer streamline or oblique-shock streamline. The 
pressure field between this streamline and the shock is assumed to be uniform. While these 
streamline locations may not be precise, the important point is that the pressure histories are 
realistic and the streamline calculations provide an accurate assessment of the nonequilibrium 
flow effects. Furthermore, for the extreme cases of the body surface streamline and the outer 
streamline, the calculations are good approximations to the results obtained from more general 
flow field calculations. 
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EQUILIBRIUM AND NONEQUILIBRIUM FLOW CALCULATIONS 

(Fig. 10) 

Calculations of the flow along selected streamlines have been made for both 
equilibrium and nonequilibrium chemical reactions. The trajectory points at which 
calculations were done are listed in the table below. Calculations have been made along 
the body surface streamline, a streamline about in the middle of the entropy layer and the 
outer streamline, or that which passes through the limiting oblique shock. The influence 
of changes in angle of attack and nose radius are also demonstrated. 

The departure from chemical equilibrium are of a different type for those 
streamlines near the body than for the outer streamlines. In the entropy layer, that is 
along those streamlines near the body surface, the flow approaches chemical equilibrium 
in the stagnation region and then departs from equilibrium in the rapid expansion of the 
flow around the nose of the body. Along the outer streamlines, the flow is very much 
like the relaxing flow behind a shock wave. The flow is chemically frozen through the 
shock wave and then approaches equilibrium over the entire body length. Thus, near the 
body surface the departure from equilibrium is a lag in the rate of recombination as the 
flow expands from the stagnation point, while along the outer streamlines there is a 
lag in the dissociation rate. 
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TRAJECTORY POINTS USED 
IN STREAMLINE CALCULATIONS 


1. 

Voo = 

4.877 km/sec, h = 

= 60.96 km 



(16 kfps) 

(200 kft) 

2. 

v OQ = 

6.096 km/sec, h = 

- 67.06 km 



(20 kfps) 

(220 kft) 

3. 

Voo = 

7.315 km/sec, h = 

* 76.20 km 



(24 kfps) 

(250 kft) 


Figure 20 



FLOW ALONG BODY SURFACE STREAMLINES 

(Fig. 11) 


Several calculations have been made for the flow along the body surface streamline, 
starting from an equilibrium stagnation point. In the altitude range of 200 to 250 kft, the 
flow along the body surface streamlines is composed mainly of N^» N, and O for both 
equilibrium and nonequilibrium flow. The principal nonequilibrium effect is the lag in the 
nitrogen atom recombination as the flow expands in the nose region. The oxygen is fully 
dissociated and remains so. The variation in temperature, density, velocity, and species 
concentrations along the streamline are illustrated in the figures below for the trajectory 
point I/p = 4. 877 km/ sec (16, 000 fps), h = 60. 96 km (200, 000 ft), and for - 1,22 m 
(4 ft), and oC = 20°. 

The temperature distribution is a direct consequence of the nonequilibrium effect 
on the N-atom recombination rate. The nonequilibrium flow temperature falls much below 
the equilibrium temperature in the nose region because the nitrogen atom recombination 
reaction becomes frozen. Then, in the fairly uniform flow on the flat aftersurface, some 
recombination takes place and the nonequilibrium temperature approaches the equilibrium 
value toward the base of the body. Notice that this approach is gradual and the flow is 
almost frozen over the entire body. The nonequilibrium effects on density and velocity 
are slight and in the direction such that there is very little effect on the mass flow rate. 



SPECIES CONCENTRATIONS, moles/g mix 


NONEQUILIBRIUM EFFECT ON SPECIES CONCENTRATIONS 





TEMPERATURE, T(°K) 


NONEQUILIBRIUM EFFECT ON TEMPERATURE ALONG BODY SURFACE 

STREAMLINE, o c= 20° 



V r n 


Figure 11(b) 
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NONEQUILIBRIUM EFFECT ON DENSITY ALONG BODY SURFACE 

STREAMLINE, c* = 20° 



V r n 


Figure 11(d) 



TRAJECTORY VARIATION OF NONEQUILIBRIUM EFFECTS ON 
THE BODY STREAMLINE FLOW 
(Fig. 12) 

Numerical solutions were obtained for the equilibrium and nonequilibrium flow 
along the body streamlines for the three trajectory points given in the table of Fig. 10. 
These calculations were done for an angle of attack of 20° and a body nose radius of 
1.22 m (4 ft). The variation of the nonequilibrium effects from one trajectory point to the 
next is one of degree rather than kind. In this series of calculations, the temperature 
distributions are all similar to that shown in Fig. 11 and so the difference between the 
equilibrium and nonequilibrium flow temperatures at the base of the body is a good 
representation of the effect of non equilibrium chemical reaction on the body streamline 
flow for these conditions. Accordingly, these temperatures are listed in the table for 
the three trajectory points at which calculations were made. The difference between the 
equilibrium and nonequilibrium temperatures is greatest at the highest altitude, highest 
velocity point because the nitrogen dissociation is greatest there and the most energy is 
frozen in chemical composition for this case. 



NONEQUILIBRIUM EFFECTS ON SURFACE STREAMLINE TEMPERATURE 


TRAJECTORY POINT 

t eq 

(°K) 

t neq 

(°K) 

1. V = 4.877 km/sec (16 kfps) 
h = 60.96 km (200 kft) 

4178 

3459 

2. V = 6.096 km/sec (20 kfps) 
h = 67.06 km (220 kft) 

4973 

3721 

3. V = 7.315 km/sec (24 kfps) 
h = 76.20 km (250 kft) 

5198 

3372 

* AT S b = 45.6 m (150 ft) 




Figure 12 


FLOW ALONG OUTER STREAMLINES 
(Fig. 13) 


The flow along the outer streamlines in the blunted-delta wing flow field is much 
like the flow behind an oblique shock wave. The reacting flow behind strong shock waves 
was extensively reviewed in Ref. 7. There several calculations were made to investigate 
scaling and distances to equilibrium. In addition, the mapping of normal shockwave flows 
to approximate wedge flows was shown to be quite accurate. The results found here are 
very much in line with what would be expected in terms of the flow behind such equivalent 
normal shocks. At the velocities and altitudes of interest the distances to equilibrium can 
be considerable. For an angle of attack of 20°, the shock angle determined by the limiting 
pressure ratio is about 2 5° and varies little from one trajectory point to the next. The flow 
behind this wave, computed here as a constant pressure streamline flow, is essentially 
frozen over the length of the orbiter. The temperature at a distance of 45. 6 m (150 ft) 
from the shock wave is shown for three trajectory points in the table below. At the lowest 
altitude of interest here, it can be seen that the flow is essentially an ideal gas flow. At 
the highest velocity and highest altitude, the difference between the equilibrium and non- 
equilibrium flow temperature is significant. Along the outer streamlines this discrepancy 
is due mainly to a lag in the oxygen dissociation rate and nitric oxide formation while the 
nitrogen remains virtually undissociated. 


NONEQUILIBRIUM EFFECTS ON OUTER STREAMLINE TEMPERATURE 

OC=20° R n = 1.22 m (4 ft) 


TRAJECTORY POINT 

a 

LU 

1- 

t neq* 


(°K) 

(°K) 

V = 4.877 km/sec (16 kfps) 
h = 60.96 km (200 kft) 

2053 

2120 

V = 6.096 km/sec (20 kfps) 
h = 67.06 km (220 kft) 

2486 

2997 

V = 7.315 km/sec (24 kfps) 
h = 76.20 km (250 kft) 

2707 

3720 


* AT s b = 45.6 m (150 ft) 


Figure 13 



NONEQUILIBRIUM EFFECTS ON THE TEMPERATURE DISTRIBUTION 

THROUGH THE SHOCK LAYER 
{Fig. 14) 

The manner in which chemical nonequilibrium effects change as attention is 
focused in turn on streamlines near the body surface and then on streamlines further out 
in the shock layer was briefly described in Fig. 10. Along the body streamline the static 
temperature of the nonequilibrium flow falls significantly below the equilibrium flow value 
as the flow expands around the nose. Along the outer streamlines, the temperature in the 
nonequilibrium flow increases toward the equilibrium value as the flow approaches chemical 
equilibrium. For one of the trajectory points, - 4.877 km/sec (16, 000 fps), h - 60.96 km 

(200,000 ft), oc = 20°, the flow along several streamlines in the shock layer has been computed. 
The results for the variation of the temperature normal to the body surface are shown in this 
figure for = 30, where is the distance along the body from the stagnation point. 

The relative values for equilibrium and nonequilibrium flow reflect the behavior of the flow 
chemistry, as described above. For this case the flow along the streamlines outside the 
entropy layer is essentially an ideal gas flow. 



NONEQUILIBRIUM EFFECT ON TEMPERATURE 
DISTRIBUTION THROUGH THE SHOCK LAYER 


ro 

h - 1 

VO 



Figure 14 
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EFFECT OF ANGLE OF ATTACK 
(Fig. 1 5 and 16) 

Streamline calculations have also been made at oc = 40°. At this angle of attack the 
flow expansion around the nose of the body is weaker. Hence there is less of a lag in the 
recombination rate for nitrogen atoms. On the afterbody the density level is higher and the 
flow approaches equilibrium much more rapidly than for oC- 20°, both along the body stream- 
lines and the outer streamlines. The temperature distribution along the body streamline and 
along the outermost streamline are illustrated below for the case Kx> = 4.877 km/sec 
(16,000 fps), h - 60.96 km (200,000 ft), and R N ~ 1.22 m (4 ft). Again these figures are 
typical of the results at the other trajectory points considered in this study. 
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NONEQUILIBRIUM EFFECT ON TEMPERATURE ALONG BODY 
SURFACE STREAMLINE, 40° 



Figure 15 
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EFFECT OF NOSE RADIUS 
(Fig. 17) 

The flow along the outer streamlines is, of course, dependent only on the distance 
from the bow wave and is not influenced by the body nose radius. On the other hand, the flow 
along the streamlines in the entropy layer are affected by changes in body scale or nose 
bluntness. In the present framework, a change in body nose radius changes the rate at 
which the pressure drops in the expansion around the blunt nose. The rate of this pressure 
drop determines the departure from equilibrium. In order to investigate the effect of nose 
radius, the body streamline calculations for oc = 20° were repeated for a body nose radius 
of . 305 m (1 ft). Again the temperature variation along the streamline has been used to 
characterize the nonequilibrium effects. The equilibrium solution and the nonequilibrium 
solutions for R.. - 0.305 m (1 ft) and R - 1.22 m (4 ft) are compared in the accompanying 
figure. As indicated by these results very little more energy is frozen in the expansion around 
the smaller nose. In addition, the rate of increase toward equilibrium is about the same in 
terms of absolute distance, the temperatures being virtually identical at a distance of 
45.7 m along the streamline. These results are typical for the range of trajectory points 
considered and nose radii in the . 305 m to 1. 22 m range. The flow along the body stream- 
line becomes frozen in the expansion around the blunt nose and recombines only a slight 
amount over the body lengths of interest. 



EFFECT OF NOSE BLUNTNESS ON BODY SURFACE 
STREAMLINE TEMPERATURE DISTRIBUTION 

= 4.877 km/sec (16 kfps) h = 60.96 km (200 kft) oc = 20° 

EQUILIBRIUM 

NONEQUILIBRIUM R N = 1.22 m 

NONEQUILIBRIUM Rm = 0.305 m 




INFLUENCE OF CHEMICAL NONEQUILIBRIUM IN THE INVISCID FLOW FIELD 

ON SURFACE HEATING RATES 
(Fig. 18) 


The calculation of the inviscid flow on the shuttle orbiter influences estimates of the 
surface heating rates through the edge conditions used in the boundary layer calculations. For 
example, consider the laminar heat transfer rate in the nose region. For a chemically frozen 
boundary layer on a noncatalytic surface, a local similarity calculation gives the heat transfer 
rate proportional to ( Ht^,~ h„) where H ^ is the total "frozen" enthalpy. Assuming that 
the body surface streamline is the edge streamline, the laminar heating rates for equilibrium 
and nonequilibrium flow in the body streamline are compared in the figure below. The results 
are shown for the case - 7.32 km/ sec (24, 000 fps), h =76.2 km (250, 000 ft), oC = 20°, 

R n = 1.22 m (4 ft). 

The assumption of frozen boundary layer flow was chosen as a point of comparison 
because the inviscid flow becomes frozen in the corner expansion and that indicates that 
the chemical state of the boundary layer will be far from equilibrium. It should also be 
noted that if the surface is completely catalytic the surface heat transfer rate depends very 
little on the chemical state of the boundary layer. Furthermore, the only way that non- 
equilibrium effects would then be important is through the transport properties. 

In order to calculate accurately the influence of the inviscid flow on the surface heating 
rate the "swallowing" of the inviscid flow by the boundary layer should be accounted for. While 
this was not done in the above calculation, it is noted that calculations of the boundary layer 
thickness indicate that the boundary layer is thin compared to the shock layer. Thus, in the 
nose region the streamlines in the entropy layer, where nonequilibrium effects are most 
severe, determine the boundary layer edge conditions. 

The entropy layer will be swallowed at some point along the body and the boundary 
layer edge conditions near the base of the body are determined by the outer streamlines of the 
inviscid flow. The ratio of ~ for nonequilibrium flow to that for equilibrium flow 
along the outer streamlines is about 1. 07 for the highest altitude, highest velocity case 
considered here. While the static temperatures are considerably different for the equilibrium 
and nonequilibrium flow along this streamline the frozen total enthalpy is not appreciably 
affected because the flow is a high Mach number flow and the static enthalpy is small 
compared to the total. Hence the nonequilibrium effects on the flow in the outer streamlines 
do not appear to have much of an effect on the surface heat transfer rate, at low angles of 
attack. 



BOUNDARY LAYER THICKNESSES, TRANSITION, AND SWALLOWING OF 

THE ENTROPY LAYER 
(Fig. 19) 

The results of the inviscid streamtube computations have been applied to estimate boundary layer 
transition, the thicknesses of the laminar and turbulent boundary layers, and the location' in the flow 
field where the entropy layer will be fully consumed by the boundary layer, i. e. entropy layer 
swallowing. Transition estimated were based on the NASA interim transition criterion. For oc = 20°, 
this criterion requires Re/S.) 10 , where all parameters are evaluated for conditions at the 

boundary layer edge. The laminar boundary layer thicknesses were estimated by using Cheng's two- 
dimensional local similarity solution,® generalized to relax the hypersonic approximations and to obtain 
a solution based on boundary layer edge conditions. For the present purposes, it was assumed that the 
boundary layer was fully turbulent immediately following transition, and the thickness estimates were 
based on the Spaulding-Chi theory9 and on procedures given by Johnson. The Spaulding-Chi theory 
was approximated by a logarithmic law, F c C f ^.o 39 [fit R e x~\~' /e > errors of less than 10%, 

and this was applied in the usual momentum integral relation for boundary layer thickness. The 
exponent in the power law velocity profile was determined from Johnson's relations, with an 
effective starting length defined to yield an approximate Blasius profile at transition. Transition and 
the boundary layer thicknesses were determined by using the various nonequilibrium streamlines in 
the entropy layer for edge conditions, and by using an equilibrium streamline processed by the 
oblique shock wave. The latter is often used in calculating boundary layer characteristics, and the 
ro present comparison provides a check on the accuracy of this procedure. 

The results of these calculations are shown in Fig. 19 for the three trajectory points considered 
here and for an angle of attack of 20°. At an altitude of 61 km (200kft), Fig. 19(aJ, the calculations show 
that the laminar boundary layer thickness is insensitive to the edge conditions, but transition is strongly 
influenced. The streamlines are labeled according to their distance from the body surface. The transition 
location varies by a factor of 2 1/2, depending on whether the nonequilibrium body streamline or the 
equilibrium oblique shock streamline is used. For these flight conditions, the latter corresponds to an 
ideal gas streamline, and there are essentially no nonequilibrium effects in that portion of the flow field. 
The thickness of the entropy layer is indicated in the figure, and it shows that the laminar boundary 
layer and transition are governed by conditions in the nonequilibrium entropy layer. The predominate 
nonequilibrium effect in this region is in the static temperature, and it differs from the equilibrium 

value by about 30%. This enters directly into the criterion for transition length through the square of 
the Mach number. 

It can be seen in Fig. 19^) that following transition, the boundary layer thickens rapidly, and it is 
thicker by a factor of 2 where the nonequilibrium entropy layer governs the flow. These boundary layer 
estimates indicate that the nonequilibrium entropy layer is fully consumed in the boundary layer at a 
distance of about 20 nose radii from the leading edge, i. e. the boundary layer development on the forward 
half of the body is governed by the nonequilibrium entropy layer and the boundary layer on only the aft 
can be estimated using conditions based on the oblique shock wave. 


ESTIMATED BOUNDARY LAYER THICKNESS 



Figure 19(a) 



(Fig. 19(b)) 


The calculated boundary layer characteristics at an altitude of 67 km (220 kft) are 
shown in Fig. 19b. These show that the thickness of the laminar layer is influenced by the 
edge conditions, and that the laminar layer develops under conditions corresponding 
approximately to the nonequilibrium body streamline. Again these results show that 
the transition length is a factor of four shorter than that predicted using equilibrium 
oblique shock conditions. The thickness of the non equilibrium entropy layer is shown 
in Fig. 19b, and it can be seen that much of the boundary layer develops under the influence 
of the entropy layer. For this flight condition, the nonequilibrium entropy layer is fully 
consumed by the boundary layer about sixteen nose radii from the leading edge. Consequently, 
nominally the forward half of the boundary layer is governed by nonequilibrium conditions in 
the entropy layer while the aft half is governed by equilibrium, or in this case, ideal gas 
conditions behind the oblique shock wave. 



ESTIMATED BOUNDARY LAYER THICKNESS 

04 = 20° h = 67 km (220 kft) V oe = 6.1 km/sec (20 kfps) R N = 1.22 m (4 ft) 







(Fig. 19(c)) 


The calculated boundary layer characteristics at an altitude of 76 km (250 kft) are 
shown in Fig. 19c. These show that the thickness of the laminar layer is influenced 
appreciably by the assumed edge conditions, differing by a factor of two for the streamlines 
considered. Again, transition is strongly and adversely effected by the nonequilibrium entropy 
layer and is predicted to occur at S^fR^ = 5 as compared with the prediction based on the 
equilibrium oblique shock conditions at S & / =17. For this flight condition, the entire 

inviscid flow field is affected by thermochemical nonequilibrium and it is not appropriate to 
calculate streamline conditions based on equilibrium flow behind the oblique shock. The 
streamline at y / = 0. 1 is processed by the oblique wave, and the nonequilibrium 
effects on boundary layer thickness can be seen by comparing the results for that streamline 
with those for the equilibrium oblique shock streamline. The entropy layer is thinner 
that at the lower altitudes and the indications are that it will be consumed by the boundary 
layer at about 12 to 14 nose radii from the leading edge. For this case, about 1/3 of 
the body length is dominated by the nonequilibrium entropy layer, and the entire inviscid 
flow field is influenced by nonequilibrium effects. 


ESTIMATED BOUNDARY LAYER THICKNESS 
oL = 20° h = 76 km (250 kft) = 7.3 km/sec (24 kfps) R N = 1 .22 m (4 ft) 



Figure 19(c) 



SUMMARY 


The degree of chemical nonequilibrium in the inviscid flow field of a blunted delta 
body typical of a shuttle orbiter was examined. This investigation was carried out within 
the framework of numerical calculations of the reacting flow along selected streamlines 
in the windward plane of symmetry. An approximate model of the pressure field was devised 
and the reacting streamline flows were computed for specified pressure histories. 

Pressure data obtained on two test shapes, a blunted delta wing and a straight-wing 
body were presented. Experimental results were presented at angles of attack of 0, 20°, 40°, 
and 60°. Streamline calculations were performed at three typical trajectory points in the 
altitude range of 60. 96 to 76.2 km (200 to 250 kft). Results were obtained for two angles of 
attack, 20° and 40°, and for two nose radii, 0.305 m and 1.22 m (1 ft and 4 ft). 

Most of the calculations were done for the 20° angle of attack and the 1.22 m nose 
radius. The results for these cases were very similar over the range of trajectory points 
considered. On the streamlines close to the body the flow approaches chemical equilibrium 
in the stagnation region and then freezes in the expansion around the spherical nose. Along 
the inner streamlines the flow is composed of essentially N^, N, and O for the conditions 
treated. Along the outer streamlines, the flow was essentially frozen for the 20° angle 
of attack. The difference between the equilibrium and nonequilibrium solutions became 
appreciable at the highest trajectory point considered. 



The effect of angle of attack was investigated by computing streamline flows for 
oC = 40° and comparing them with the oL = 20° results. The increase in density level on 
the windward surface was sufficient to promote an approach to chemical equilibrium along 
both the inner and outer streamlines in the shock layer. The calculations carried out for 
R = 0.305 m (1 ft) indicated very little effect of nose radius or scale on the degree of 
nonequilibrium. This is because the flow along the body streamline freezes even for larger 
values of the nose radius. The streamline outside the entropy layer are independent of 
nose bluntness effects. 

The edge conditions of the laminar boundary layer in the nose region are determined 
by streamlines which pass through the nearly normal portion of the bow shock. On the other 
hand the streamlines which pass through the limiting oblique shock determine the edge 
conditions of the boundary layer toward the base of the body. The calculations done here 
for an angle of attack of 20° show a fairly large fraction of the total enthalpy, e.g. about 
30%, can be frozen in the chemistry along the body streamlines. On the other hand, in the 
flow along the outer streamlines the chemistry is nearly frozen but accounts for less than 
10% of the total enthalpy. As the angle of attack increases the fraction of the total enthalpy 
accounted for by dissociation increases, but so does the rate of dissociation. At an angle 
of attack of 40°, the flow in the outer streamlines is near equilibrium over most of the 
length of the vehicle. 



The influence of nonequilibrium effects on the boundary layer has been estimated 
using approximate boundary layer calculations. These show that for all three trajectory 
points, the laminar boundary layer development is governed by the nonequilibrium 
entropy layer. In each case, transition is determined by conditions in the nonequilibrium 
entropy layer and not by conditions behind the oblique wave. In each case the effect was 
detrimental in that transition occurred substantially earlier than predicted on the basis 
of equilibrium flow behind the oblique shock. In all cases covered, the nonequilibrium 
entropy layer governs the boundary layer development on the forward 1/3 to 1/2 of the 
body. The boundary layer development on the aft portions is governed by equilibrium flow 
behind the oblique shock at the two lower altitudes, and by nonequilibrium flow behind 
the oblique shock wave at the highest altitude. 

There are several features of the nonequilibrium flow over a shuttle orbiter which 
deserve mention before concluding. First, it should be noted that the flow in the windward 
plane of symmetry does not necessarily represent the extreme so far as chemical non- 
equilibrium effects are concerned. Along streamlines which enter the nearly normal 
portion of the bow shock on the windward side but are then swept around the edge of the 
body to the leeward side, the flow undergoes a much greater expansion than in the plane 
of symmetry. Consequently, the difference between the equilibrium and nonequilibrium flow 




properties and composition along such streamlines could be greater than along those treated in this 
study . 

Although the presence of ionization was allowed for in the calculations the results for the elec- 
tron concentration have not been presented in detail. However, the effect of nonequilibrium in the 
ionization and deionization reactions is even more pronounced than on the neutral chemistry. For the 
range of conditions considered here the electron number densities were between 10^ cc ^ and 10^ cc ^ 
in the nose region and between 10^ cc ^ and 10^ cc ^ at the base of the body. 

The results of the present study indicate that nonequilibrium effects will persist at altitudes 
below 60.96 km (200,000 ft) for those streamlines near the body surface. Also nonequilibrium effects 
along the outer streamlines could be important at higher angles of attack at altitudes above 76.2 km 
(250,000 ft). 

For the range of trajectory points considered here the nonequilibrium effects on a blunted delta 
shape at an angle of attack of 20° are much like those in the flow on blunted slender cones. Non- 
equilibrium effects are most important in the entropy layer. However, significant effects on the tem- 

/’ 

perature and composition exist throughout the flow field in portions of the orbiter entry trajectory. 
The influence of the nonequilibrium effects on the boundary-layer edge conditions have been discussed. 
However, the nonequilibrium effects on the viscous portions of the flow field have not been treated 
beyond this . Furthermore the results of the present work indicate that the flow in the boundary layer 
will be far from equilibrium in the altitude range above 60.96 km (200,000 ft). 
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